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1. Introduction 

The asymptotic freedom of QCD guarantees a small coupling constant g at large temperatures 
T. While observables can be expressed in a generalized power series in g, the loop expansion is 
not applicable to an arbitrary order in g, because of the so-called "infrared wall", as pointed out by 
Linde [[Tj] (see also ref. [§]). For every observable there exists an order of the perturbative expansion 
to which an infinite number of Feynman diagrams contributes. For the pressure this order is g 6 T 4 . 

No way of resumming these infinitely many diagrams has been found, so a different approach 
is needed. The problem arises due to infrared divergences in the dynamics of zero Matsubara 
frequency modes of gauge fields. Because these modes are three-dimensional (3d) we can construct 
an effective 3d pure gauge theory called Magnetostatic QCD (MQCD) which accounts for the 
non-perturbative contribution [Q, f|, §]. QCD and MQCD can be matched to each other by using 
perturbation theory. 

To obtain the non-perturbative contribution we perform lattice measurements in MQCD [§]. 
The observable we consider is the elementary plaquette expectation value. The theory being super- 
renormalisable, one can match the lattice regularization scheme exactly to the MS scheme. This 
requires a perturbative 4-loop computation on the lattice, however, which has not been completed 
yet: the missing ingredient is specified below. (A certain perturbative 4-loop computation in full 
QCD remains also to be carried out.) 



2. Relation between MS and lattice regularization schemes 

The euclidean pure SU(/V C ) Yang-Mills action reads 

S E = f d d X Jf E , jgfe = -L£Tr[/^], 

■' 2 <?3 k,l 



(2.1) 



iA k ,A k =A a k T a , 



where d = 3 — 2s, gj is the gauge coupling, k,l = 1, . . . ,d, Fu = i[Dk,Di],Dk = <9& 
and T a are hermitean generators of SU(Af c ) normalised such that Tr[T a T b ] = 8 ab /2. The vacuum 
energy density is defined as 



/ms — 
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J@A k exp(-S E ) 



(2.2) 
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where V denotes the (i-dimensional volume. The use of the MS dimensional regularization scheme 
removes any l/e poles from the expression. In fact, using dimensional regularization the perturba- 
tive result vanishes, because there are no mass scales in the propagators. However, for dimensional 
reasons, the non-perturbative form of the answer is 
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where d& = N% — 1. The logarithmic term has been calculated by introducing a mass scale for 
gluon and ghost propagators and sending — > after the computation [[7|, . 

Using standard Wilson discretization, we can write the same theory on the lattice as 
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(2.4) 
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where Py is the plaquette, a is the lattice spacing and j3 = 2N C / (ag 2 ). Hence the continuum limit is 
taken by j8 — > °°. Dimensionally, the vacuum energy density consists of terms of the form g 2n a n ~ 3 . 
Thus, approaching the continuum limit, we can relate f a and as follows: 



Af = f a -fus (2-5) 

ay 



C i 3 ( ln + C 'i V C 2 4 + C 3 - + C 4§3 f l n i + + 0{&a) . (2.6) 

a V a si J a a \ fl/i 7 



Taking derivatives of eq. (^3|) with respect to g\ and using 3d rotational and translational 
symmetries on the lattice, we obtain the master relation 



</ ^G=limy3 4 ((l--iTr[P]) a 



(4tt) 4 u p_»=^ 1/ N, 



c\ c 2 c 3 c 4 , 
]8 + F F 4) 



(2.7) 



This quantity may be called the finite part of the gluon condensate in lattice regularization (in 
certain units). 

The values of the constants c\ , . . . , c* 4 are trivially related to those of C\, . . . ,C' A . When N c = 3, 
the numerical values for the c ; 's are 

Cl = -A a 2.666666667, (2.8) 
c 2 = 1.951315(2), (2.9) 
c 3 = 6.8612(2). (2.10) 

Here ci results from a straightforward 1-loop computation, while C2 and C3 have been calculated 
in refs. [^] and [10], respectively. Because there is no p-dependence in /„, the value of C4 is 
determined by in eq. ( |23| ). Consequently, 

c 4 « 2.92942132. (2.11) 

Because the constant c\ = C\ — 1/3 — 21n(2A/c) is still unknown, we are not able to fully 
determine Bq. We can, however, determine the non-perturbative input needed for it. In order 
to evaluate c' A a 4-loop lattice perturbation theory calculation is required; alternatively, it can be 
evaluated non-diagrammatically by means of numerical stochastic perturbation theory 



3. Lattice measurements 

We need the plaquette expectation value (1 — ^Tr[P]) a as a function of j8 so that the extrapola- 



tion j3 — > 00 in eq. (2.7) can be carried out. For each j8 the infinite- volume extrapolation is needed. 



Due to the non-perturbative mass gap of the theory, the finite-volume effects are exponentially 
small if the size of the box L = Na is large compared with the inverse confinement scale g^ 2 . In 
practice finite- volume effects are invisible as soon as f5/N < 1. This is demonstrated in Fig. |l[ In 
Fig. H| the infinite- volume extrapolated values of (1 — ^Tr[P]) a are plotted as a function of 1/jS. 
A major difficulty in the simulations is the significance loss caused by the subtractions in 



eq. (2.7). The dominant term C[/f5 is about six orders of magnitude larger than the effect we are 
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Figure 1: Finite-volume values for /3 4 {(1 - jTr[P]) a - [ci/j3 + c 2 //3 2 + c 3 //3 3 +c 4 rnj3/j3 4 ]} as a function 
of the box size. The leftmost symbols indicate the infinite-volume estimates, obtained by fitting a constant to 
data in the range j3 /N < 1 (blue diamonds). Points denoted by red circles are omitted from the extrapolation. 
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Figure 2: The plaquette expectation value plaq = (1 — |Tr[P]) u as a function of 1//3. Statistical errors are 
much smaller than the symbol sizes. The solid line contains the four known terms ci//3+C2/j3 2 + C3/j3 3 + 
C4 In [i /j3 4 . The effect we are looking for is the difference between the data and the line. 



interested in, namely ~ 1//3 4 , if j8 ~ 100. Therefore the relative error of our lattice measurements 
should be smaller than one part in a million. The effect of the subtractions is illustrated in Fig. ||. 

Given the infinite-volume limits, we extrapolate the data to the continuum limit, j3 — > °°. In 
Fig. | we show two functions: /3 4 {(1 - jTr[P]) Q - [ci/j3 + c 2 /j8 2 +c 3 /j3 3 ]} and /3 4 {(1 - ^Tr[P}) a - 
[ci//3+C2//3 2 + C3/j8 3 + C4lnj8//3 4 ]}. Even the 4-loop logarithmic term is visible in the data. For 
1/P < 0.01 the significance loss grows rapidly and the error bars become quite large, so that these 
data points have little effect on the fit. 
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Figure 3: The significance loss due to the subtraction of divergent lattice contributions. Here again plaq = 
(1 — ^Tr[f]}„, and the symbols c, indicate which subtractions of eq. (2.7) have been taken into account. 




Figure 4: The infinite-volume extrapolated data. The effect of the 4-loop logarithmic divergence is to cause 
additional upwards "curvature" in the upper data set. The solid blue line gives the continuum extrapolation. 
Points with lighter colors have so large errors that they are insignificant as far as the fit is concerned. 
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The continuum extrapolation is carried out by fitting a function d\ + ^/jS + d^/fi 2 to the 
infinite-volume extrapolated data, from which all the divergences ({ci ,C2,C3,C4}) have been sub- 
tracted, in the range 0.01 < 1/j8 < 0.10. The fitted values are d\ = 20.0(7), d 2 = 86(24) and 
^3 = 909(192) with # 2 /dof= 5.8/6. The error limits are the projections of the 68% confidence 
level contour onto the various axes. The systematic errors from the effect of higher order terms are 
inside these errors. 



Substituting this to eq. (2.7) we obtain the final result, 



/43 157 „\ 47T 4 
Bg + ( 15 " 768 71 ) C ' 4 = ~W X 20 -° (7) = 10 - 7(4) - (3 " 1} 



4. Conclusions 



We have studied the expectation value of the elementary plaquette in 3d pure SU(3) theory and 
outlined how the 3d vacuum energy density in the MS scheme can be extracted from it. However, 



to achieve this, the constant c' 4 should be determined (cf. eqs. (|2J|), (3.1)). This can be accom- 
plished by a 4-loop matching computation, with the techniques discussed in ref. [11]. The full 
QCD pressure of order g 6 T 4 can be obtained by a further 4-loop matching computation. 
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